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Recently, it was shown that giant vortices with arbitrarily large quantum numbers can possibly 
be created in dilute Bose-Einstein condensates by cyclically pumping vorticity into the condensate. 
However, multiply quantized vortices are typically dynamically unstable in harmonically trapped 
nonrotated condensates, which poses a serious challenge to the vortex pump procedure. In this 
theoretical study, we investigate how the giant vortices can be stabilized by the application of 
a Gaussian potential peak along the vortex core. We find that achieving dynamical stability is 
feasible up to high quantum numbers. To demonstrate the efficiency of the stabilization method, 
we simulate the adiabatic creation of an unsplit 20-quantum vortex with the vortex pump. 

PACS numbers: 03.75.Lm, 03.75.Mn, 67.85.Fg 



I. INTRODUCTION 

Bose-Einstein condensation in trapped gases of alkali- 
metal atoms was observed experimentally in 1995 (THl]. 
A few years later, those pioneering experiments were fol- 
lowed by the creation of singly quantized vortices 0, H[ 
and vortex lattices 041] in such systems. Since quantized 
vortices manifest the long-range phase coherence of the 
condensates and are inherently connected with superflu- 
idity, there has been wide interest in studying their prop- 
erties in trapped Bose-Einstein condensates (BECs) [T(| ■ 
In particular, the stabilit y o f vortices has been the sub- 
ject of intensive research fllH23| . 

In principle, a quantized vortex in a BEC can carry 
any number of circulation quanta. However, it is well 
known that typically a vortex with a winding quantum 
number greater than unity has a higher energy than the 
corresponding number of separated single-quantum vor- 
tices, which means that multiquantum vortices have a 
propensity to split into single-quantum vortices [2^ - l3l| . 
Although this result is generally true only in an infinite 
homogeneous system, it still holds in finite-sized BECs 
for a majority of trap geometries and particle numbers. 
Being able to create vortices with large winding numbers 
would provide access to novel vortex-splitting patterns 
beyond the typical linear chain that prevails in the decay 
of two- and four-quantum vortices [2{| [28|. Because of 
the distinct nature of the different splitting patterns pre- 
dicted for highly quantized vortices [30| . observing the 
decay of such states would enable a lucid comparison 
between theory and experiment. It has also been spec- 
ulated that the splitting of multiquantum vortices may 
create necessary conditions for the initialization of super- 
fluid turbulence [HI]. 

The energetic instability renders the creation of mul- 
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tiquantum vortices challenging but not impossible. In 
fact, energetically unstable states can be quite long-lived 
since the relaxation to lower-energy states necessitates 
dissipation which comes mainly in the form of noncon- 
densed atoms. Temperatures in typical experiments are 
low enough such that the noncondensed component of 
the gas is negligible and, consequently, relaxation is slow. 
Therefore, methods that do not rely on the relaxation of 
condensate dynamics are especially well-suited for pro- 
ducing multiquantum vortices. In topological phase im- 
printing |33l436| . the condensate remains in its instan- 
taneous eigenstate, and the process is insensitive to the 
exact rate of change and path of the control fields. 

Although topological phase imprinting was the first ex- 
perimentally realized means to create multiquantum vor- 
tices 37 1, it was only recently [38j when it was shown to 
be applicable in creating giant vortices with essentially 
arbitrarily large winding numbers. In this vortex pump, 
a fixed amount of vorticity is added to the condensate 
in each adiabatic pumping cycle, and thus in principle, 
arbitrarily large vorticities can be reached. However, a 
serious challenge is posed by the dynamical instability of 
the giant vortices that becomes more pronounced as the 
vorticity increases 0, [23[ . Dynamical instabilities can 
lead to dissociation of the vortex even in the absence of 
dissipation [2^ - [3l| . and therefore their effect cannot be 
disposed of by reducing temperature. 

Motivated by the above considerations, we study how 
vortices with large winding numbers can be made dynam- 
ically stable in nonrotated harmonically trapped BECs. 
As our method of choice, we investigate the effect of ap- 
plying a repulsive Gaussian-shaped plug potential along 
the symmetry axis of the trap. In the vortex pump, the 
plug not only serves to stabilize the vortex but also pre- 
vents unwanted spin flips in the central region of the trap. 
The plug potential can be realized in experiments by a 
tightly focused far-blue detuned laser beam [|| |39| - 41 
as was done already in the pioneering work of Davis et 
al. 1. 
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This article is structured as follows. In Sec. |TT1 we 
present the mean-field theory of BECs and relate the 
computational parameters of the plug potential to the 
properties of the laser. We also describe how the vortex 
pump is modeled. In Sec. IIII1 we present our numer- 
ical results and show that vortices up to high winding 
numbers can be stabilized with experimentally achiev- 
able plug potentials. We demonstrate the efficiency of the 
plug by simulating the creation of an unsplit 20-quantum 
vortex with the vortex pump. Section IIVI is devoted to 
discussion. 



Small-amplitude oscillations about the statio nary 
states play an important role in the study of BECs [4J] . 
We decompose the order parameter as 

*(r, 4>, t) - e-^ h [*(r, 0) + X (r, 0, t)\ , (4) 

where we assume the oscillatory part x to have a small 
L 2 norm compared with vp. In the Bogoliubov theory, 
one seeks a solution in the form 
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II. THEORETICAL AND NUMERICAL 
METHODS 

We consider a BEC in a nonrotated, cylindrically sym- 
metric harmonic trap and restrict the analysis to the 
zero-temperature limit, thereby ignoring the effects of the 
thermally excited atoms [42| . Experiments with dilute 
BECs can usually be carried out at temperatures where 
this approximation is justified (43j . In the stability anal- 
ysis below, we assume that the spin degree of freedom 
of the BEC is fixed by the Zeeman coupling to a strong 
external magnetic field. Under these circumstances, the 
scalar order parameter field ^ of the BEC is described 
by the time-dependent Gross-Pitaevskii (CP) equation, 

ihd t ^{r,cp 1 t)= [H + g\y(r^,t)\ 2 ] *(r,<f>,t), (I) 

where % denotes the single-particle Hamiltonian and the 
atom-atom interaction strength g is related to the vac- 
uum s-wave scattering length a, the atomic mass m, and 
the axial harmonic oscillator length a z = sjKj (muj z ) by 
g = \/8nh 2 a/ (ma z ). In Eq. (TTJ), we have assumed that 
the condensate is pancake-shaped, i.e., that the harmonic 
trapping frequency in the axial direction is much greater 
than in the radial direction, lj z 3> 0J r , which has enabled 
us to factor out the z dependence of the full order pa- 
rameter as \&f u ii(r, t) = ^{r, 0, t) exp[— z 2 / (2a 2 )]/ \J~Kd 2 . 
The single-particle Hamiltonian is defined by 
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~dr + -A 



V(r), 
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where the potential function includes the possible optical 
plug potential, V(r) = mco 2 r 2 /2 + Vpi ug (r). The order 
parameter is normalized such that J |\&| r drd0 = N, 
where N is the number of condensed atoms. 

Stationary states of the system satisfy the time- 
independent GP equation, which is obtained from Eq. (fl]) 
with the replacement ihd t — > /i, where /i is the chem- 
ical potential. Here, the stationary states are chosen to 
be axisymmetric vortex states with a winding number k, 
implying that the order parameter can be written in the 
form 



where n(r) is the areal particle density. 



(3) 



where the complex- valued functions u q and v q are the 
quasiparticle amplitudes corresponding to the index q. 
Each quasiparticle mode is also characterized by an in- 
teger l q that determines the angular momentum of the 
mode with respect to the condensate. By substituting 
Eqs. ((H) and ([5]) into Eq. (fTJ), neglecting terms super- 
linear in x, and noting that ^{r, </>) satisfies the time- 
independent GP equation, we arrive at the Bogoliubov 
equations 
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where 



2m 



Or 



+ V(r)-n + 2gn(r). (7) 



The Bogoliubov excitation spectrum {uj q } can be used 
to classify the stability of the corresponding stationary 
state. If the spectrum contains excitations with a positive 



norm 



/[It 



r dr but a negative eigenfrequency 



co q , the stationary state is energetically unstable. States 
that support modes with nonreal eigcnfrcqucncies are re- 
ferred to as dynamically unstable, because the amplitude 
of a small perturbation associated with the excitation 
of a complex-frequency mode initially evolves exponen- 
tially in time [Eq. ((SJ)]. In the case of a multiquantum 
vortex, dynamical instability typically signifies that the 
vortex is unstable against splitting into singly quantized 
vortices (H^n]]. 

In the numerics, we measure length in units of the 
radial harmonic oscillator length a r = \J hj {muj r ) and 
energy in units of hw T and normalize the dimensionlcss 
order parameter to unity. With these scalings, the dimen- 
sionless interaction strength becomes g = \Z8nNa/a z . In 
experiments, g lies typically between 10 2 and 10 6 . 



A. Vortex stabilization with an optical plug 

As our first topic, we study the application of a re- 
pulsive plug potential as a means of dynamically stabi- 
lizing vortices with large winding numbers. In the case 
of giant vortices, dynamical instabilities correspond to 
deformations of the vortex core, and thus the quasipar- 
ticle amplitudes of dynamically unstable excitations are 
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nonvanishing within the core region. Consequently, the 
instabilities of an axisymmetric giant vortex can be sup- 
pressed with a sufficiently strong and wide plug potential 
along the symmetry axis of the trap, since it increases the 
energy of such excitations. We assume that the plug has 
a Gaussian profile, i.e., 



Vpiug(r) = Ae 



1 I ,2 
jd 



(8) 



where A denotes the amplitude and d is the beamwidth. 
The dynamical stabilization of a giant vortex was briefly 
investigated in Ref. [30| in the case of an infinitely hard 
step-function potential, and the width of the step re- 
quired to stabilize the vortex was found to be somewhat 
smaller than the size of the vortex core. 

In order to relate the parameters in Eq. (|8|) to ex- 
perimental quantities, let us derive the expression from 
atomic properties. If the laser is operated at a fre- 
quency u> and the resulting electric field is denoted by 
E, the potential experienced by an atom is given by 
Vpiug = — a(w)E 2 /2, where a(uj) is the dynamical po- 
larizability and (•) denotes the time average. The spatial 
dependence in Eq. © follows from the spatial profile 
of the electric field, which for a focused laser beam is 
taken to be Gaussian. The polarizability is given by the 
Kramers-Heisenberg formula [44j 
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where e is a unit vector in the direction of the electric 
field, d is the electric dipolc-moment operator, and we 
label the ground state by g and the excited states by e. 
The transition energies are expressed as huj eg — e e — e g > 
0. We assume that the light is detuned far from the 
relevant atomic resonance at ujq such that the detuning 
frequency ui — wo is much larger than the natural decay 
rate of the corresponding excited state. For the two most 
common BEC species, 23 Na and 87 Rb, ujq corresponds to 
the n 2 Si/2 — > n 2 P 3 /2 transition, i.e., to the D2 line, with 
angular freq uencies given respectively by 3.20 x 10 15 Hz 
(n = 3) [Htl and 2.41 x 10 15 Hz (n = 5) @. Thus, we 
limit the summation in Eq. ([9]) to the single hyperfine 
manifold of the n 2 P 3 / 2 level and neglect the hyperfine 
splitting between the different transitions, which yields 
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where we have labeled the remaining sum in terms of 
the effective dipolc moment d c ff. The actual value of 
|d e fj| depends on the atomic ground state and the po- 
larization of the light but can, nevertheless, be readily 
evaluated [U HH]. In the case of 87 Rb atoms in an 
F = 1 hyperfine state and 7r-polarized light, one obtains 
|d e fj| = 2.44eoaB, where e is the electron charge and ae 
is the Bohr radius. If we further assume that the laser 
operates at 660 nm (a typical experimental value) and 



relate the maximum electric field to the power P of the 
laser, |E| 2 lax = 2P/ (wceod 2 ), we find 



VpiugM = fc B x 73 /jlK 
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To study the dynamical stabilization of giant vortices 
numerically, we solve the stationary GP equation and 
the Bogoliubov equations for different values of the wind- 
ing number k, interaction strength g, plug amplitude A, 
and beamwidth d and assess the dynamical stability of 
the corresponding stationary state. We use finite dif- 
ference methods and solve the stationary GP equation 
using successive over-relaxation. The Bogoliubov equa- 
tions are solved using the LAPACK numerical library in 
MATLAB lil. 



B. Optical plug in vortex pumping 

As our second topic, we demonstrate that by utiliz- 
ing a sufficiently strong plug potential, the vortex pump 
method [38| can be used to create an unsplit giant vor- 
tex with a very large winding number. To this end, we 
study the temporal evolution of a BEC during vortex 
pumping in a case where the harmonic trap is combined 
with a strong Gaussian potential of the form of Eq. . 
Since the vortex pump makes explicit use of the spin de- 
gree of freedom of the BEC, we no longer assume the 
atomic spins to be fixed. Instead, we consider a spinor 
BEC with a hyperfine spin F = 1 and model it with the 
time-dependent spin-1 GP equation [50L l5l| 



ihdt^ir, <j>,t) 



H + m9FB(t) ■ T + c &$ 
-c 2 (&F$) .^(r.&t), (12) 



where T = {F x ,F yi F z ) 
tors of the spin rotation 



consists of the genera- 
group SO (3) and * = 



(^i, ^0, ^-i) 1 is the three-component spinor order pa- 
rameter written in the eigenbasis of T z . Moreover, /xb 
is the Bohr magneton, gp is the Landc g factor, and 
Co = ^/Snfi 2 (clq + 2a 2 ) / (3ma z ) and c 2 = V8nh 2 (a2 — 
ao)/ (3ma z ) are the coupling constants related to s-wave 
scattering lengths ag and a 2 for different spin channels. 
As in Eq. ([1]), we have again assumed that the BEC is 
pancake-shaped and integrated out the z dependence of 
the order parameter. 

In the vortex pump, the spin degree of freedom of the 
condensate is controlled locally by slowly tuning the ex- 
ternal magnetic field B(t) in a cyclic manner such that 
the system acquires a fixed amount of vorticity in each 
control cycle. The efficient operation of the pump re- 
quires that sudden spin flips due to Landau-Zener tran- 
sitions are insignificant, and therefore, the control cy- 
cle should be sufficiently adiabatic [52| . In our simula- 
tion, we use a cycle which is identical to the one em- 
ployed in Ref. [381 ] . consisting of a homogeneous bias 
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Figure 1: (Color online) (a) Control cycle of the vortex 
pump simulation in the magnetic-field parameter space. Here, 



B' 



, g \/r and B' h 



\Bh\/r denote the radial derivatives of 



the quadrupole and hexapole field magnitudes, (b) Temporal 
changes of the control parameters during the pumping cycle 
given in Eq. (|13[) . The common maximum value of B' q and 
B' h is denoted by B®. 



field in the z direction and alternating quadrupole and 
hexapole fields in the xy plane. Denoting the quadrupole 
and hexapole fields by B° = B®r [cos(0)x — sin(0)y] and 
B° = £?J?r[cos(20)x — sin(2</>)y] , the exact form of the 
control cycle becomes 
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Figure 2: Maximum imaginary part of the eigenfrequencies for 
a fO-quantum vortex as a function of the interaction strength 
g for different values of the plug amplitude A given in the 
inset. The width of the plug is set to d = 3a r . The imaginary 
parts vanish completely for A > 26.0ftu r - The dimensionless 
interaction strength is given by g = \^8TrNa/a z , where N 
is the particle nu mber, a the s-wave scattering length, and 
a z = ^Jhj (mujz) the axial harmonic oscillator length. 



III. RESULTS 
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< x < 1, is such that the spin at a distance of 5a r from 
the z axis is reversed with a constant speed. Here, T\ is 
the ramping time of the multipole fields, T 2 is the bias- 
field inversion time, and T = 4T\ + 2T 2 is the total period 
of the cycle. The control cycle is visualized in Fig. [TJ In 
order to make the effect of the plug potential transparent, 
we set the system parameters identical to those used in 
the original simulation of Ref. 38| except for the plug 



potential, which we choose to be significantly stronger. 
Thus, we aim at achieving a significantly higher winding 
number than the value k = 8 reported in Ref. [38j . Before 
beginning the pumping, we use the method of successive 
over-relaxation to find the ground state of the condensate 
with the magnetic field set to its initial configuration, 
B(i = 0) = —B z z. Equation (|T2"j) is then numerically 
integrated using the Strang splitting scheme. 



A. Vortex stabilization with an optical plug 

By solving the stationary GP and Bogoliubov equa- 
tions numerically, we have determined the values of the 
plug amplitude A and width d that are sufficient to ren- 
der the multiquantum vortex states dynamically stable 
at different winding numbers k and dimensionless inter- 
action strengths g = \/8irNa/a z . The value of g has been 
restricted to the range < g < 2000, but this should not 
limit the generality of our results. 

Figure [2] shows the strongest dynamical instability of 
a 10-quantum vortex as a function of g, < g < 1000, 
for different values of the amplitude A. In the absence of 
the plug, the vortex is dynamically unstable in the whole 
interval < g < 1000, but as the strength of the plug 
increases, regions of dynamical stability begin to appear, 
until eventually the unstable modes vanish completely. 

In Fig. [31 we plot, as a function of k, the limiting plug 
amplitude A m i n above which the K-quantum vortex is 
dynamically stable in the whole interval < g < 2000. 
Here, the beamwidth is fixed at d = 3a r , but quanti- 
tatively similar behavior is found for other values of d 
as well. The stabilizing amplitude A m \ n increases expo- 
nentially with k for sufficiently large winding numbers, 
and ordinary least-squares fitting yields the expression 
= 13ftw r exp(0. 072k). Roughly speaking, the expo- 
nential dependence follows from requiring that the plug 
potential exceeds a certain limiting strength Vq inside 
the whole core region r < r c , where r c denotes the vor- 
tex core radius. The plug is weakest at r ~ r c , and thus 
the limiting condition becomes ^4 m i n exp(— r^/d 2 ) = Vq. 
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Figure 3: (Color online) Minimum amplitude A m i n [see 
Eq. (jlll) ] required to stabilize vortices with a given winding 
number n. The black dots indicate values that stabilize the 
vortices in the whole interval < g < 2000, and the dashed 
blue line is the least-squares fit j4 m i n = 13fko r exp(0.072fi;). 
The green diamonds show the values that stabilize the vortex 
for fixed g = 250, and the dotted red line represents the fit 
Anin = 6.2hu r exp(0.065fc). In both cases, the width of the 
plug is set to d — 3a r . 



3.5 
3 

2.5 
2 



a 

£ 1.5 



t3 



1 

0.5 









10 20 30 40 50 



Figure 4: (Color online) Minimum beamwidth rf m j n [see 
Eq. (flip) ] required to stabilize vortices with a given winding 
number k in the whole interval < g < 2000. The dot- 
ted curve is a two-parameter fit to a square-root function for 
k > 10, given by d m i u /a r — 0.52y/ti — 6.8. The amplitude of 
the plug is set to A = 200huj r . 



It was found in Ref. [23[ that for a sufficiently large n, 
r c /a r oc t/k. Neglecting the dependence of Vq on k, 
we get log^min cx k [53[. The irregular behavior at 
k < 10 is explained by noting that the dependence of 
r c on k deviates significantly from the square-root form 

if ,g>K 2 m. 

The obtained values of ^4 m i n are small enough to be 
realizable with commercial lasers: If we use the values 



25 kg corresponding to 
= 2tt x (8,90) Hz EL 



a = 4.7 nm and m = 1.44 x 10" 
87 Rb atoms choose (uj r ,uj z ) 

and use the plug potential in Eq. (fTTj) . we find that d 
3a r ss 11 ^m and that ^4 m i n = 454ftu; r at k = 50 is 
obtained with the laser power P = 0.31 mW. Values of 
this order have been used in various experiments [H, 0, 

MM- 

The smallest beamwidth c? m i n required to stabilize the 
multiquantum vortices with < g < 2000 and a given 
winding number k is displayed in Fig. |4] as a function of 
k. The amplitude is fixed at A = 200fi-w r . The minimum 
width d m i n increases as a square-root function of k. This 
can again be understood by assuming that the plug is 
stabilizing when its strength exceeds the threshold value 



Vo for r < r c . Thus, we require Aexp(—r^/df nin ) = Vq, 
which yields d mm oc r c oc ^/k. 



B. Optical plug in vortex pumping 

To demonstrate the efficiency of the stabilizing plug, 
we have computed the time evolution of the spinor or- 
der parameter ^ from Eq. (|12[) during vortex pumping. 
The optical plug has been chosen to have the amplitude 
A = 100fio; r and width d = 3a r , and thus it is consider- 
ably stronger than that used in Ref. [38| (A = 10fkj r 
and d = 2a r ). Otherwise, the system is identical to 
the one considered in Ref. |38|, where unsplit vortices 
up to the winding number k = 8 were reached. More 
specifically, the dimensionless interaction parameters are 
c = mNco/h 2 = 250 and c 2 = mNc 2 /h 2 = -0.01c 
(which corresponds to 87 Rb), the magnetic field strengths 
read = 40hui r / (^b\9f\) and B® = hw r j (a r ^B\gF\), 
and the multipolc-field-ramping and bias-field-inversion 
times are given by T± = 10/uj r and T 2 = 160/cj r such 
that the cycle period is T = 360/o; r [see Eq. (fl3"|)]. 
The Lande g factor gp is taken to be negative. As- 
suming again that the BEC consists of 87 Rb atoms and 
(av,w z ) = 27rx (8, 90) Hz [4l|, the simulation parameters 
correspond to a r ps 3.8 fj,m, iV k 10 4 , ~ 46 nT, and 
B° w 0.3 nT/^m. 

In Fig. [5l we present the squared modulus and complex 
phase of the order parameter component at integer 
multiples of the pumping period T. The accumulation 
of two quanta of vorticity per cycle is clearly visible in 
the phase field. Figure [6] shows the corresponding time- 
dependence of the total axial angular momentum (L z ) of 
the BEC. The notable deviation of (L z ) from the ideal 
value 2kNh after k rj 8 pumping cycles is due to excita- 
tions away from the instantaneous eigenstate the pump is 
operated in. These excitations are imperceptible in Fig. [5] 
during the nine first cycles, so that a clear 18-quantum 
vortex state is observed at t — 9T. In fact, a symmetric 
20-quantum vortex is obtained in the middle of the ninth 
cycle at t = 8.5T. We have also checked that the plug 
amplitude can be ramped down without destroying the 
giant vortex. 

After t = 9T, the dissociation of the giant vortex be- 
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Figure 5: (Color online) Areal particle density (upper panels) 
and complex phase (lower panels) of the component of 
the spinor order parameter in the xy plane at integer multi- 
ples of the vortex pumping period T — 360/oj r , with a plug 
potential of width d = 3a r and amplitude A = 100hu} r . The 
field of view in each panel is 18 a r x 18 a r . 

gins in spite of the plug, causing single vortices to move 
out of the condensate. This is also manifested in Fig. [5] 
by the decrease of the axial angular momentum. More- 
over, by using initial states with higher vorticity, we have 
verified that the splitting after t = 9T is not caused by 
the accumulation of numerical errors in the temporal evo- 
lution and that a higher winding number is reached by 
further increasing the strength of the optical plug. We 
have also observed that the abrupt removal of the plug 
and the addition of weak random noise to the state at 
t = 8.5T or t = 9T eventually results in the splitting of 
the 18- or 20-quantum vortex with a fourfold-symmetric 
splitting pattern introduced in Ref. [30| . 

IV. DISCUSSION 

In this article, we have studied the dynamical stabi- 
lization of giant vortices in harmonically trapped BECs 
by applying a Gaussian-shaped repulsive plug potential 




Figure 6: (Color online) Axial angular momentum of the con- 
densate as a function of time during the vortex pump sim- 
ulation presented in Fig. [5] The inset shows the particle 
density and phase of the 20-quantum vortex (field of view 
18 a r x 18 a r ) at the time instant t = 8.5T marked with the 
red cross. 



along the symmetry axis of the trap. We found that vor- 
tices with large winding numbers can be stabilized with 
plug profiles that should be routinely achievable with 
commercial lasers. Although the detailed behavior of the 
dynamical instabilities as a function of the plug param- 
eters turned out to be complicated (Fig. [2]) , the overall 
criterion for the stabilization could be explained in sim- 
ple terms of filling the vortex core with a sufficiently high 
potential barrier. 

We also performed a simulation of the vortex pump 
which indicated that a giant vortex can be created by 
pumping if a sufficiently strong optical plug is utilized 
and the temperature is kept low enough such that dissi- 
pation effects due to the thermal gas are negligible. In 
the simulation, the duration of the control cycle was kept 
constant throughout the creation process. In fact, it is 
possible to gradually increase the pumping speed as vor- 
ticity accumulates into the BEC and still retain the adia- 
baticity of the process [23[ . By making use of this possi- 
bility, one could increase the winding number beyond the 
value k — 20 reached in our simulation before the vortex 
splits. Furthermore, one could also employ a different 
pumping scheme (lil . [55| which increases the vorticity of 
the F = 1 condensate by 4h/m per cycle instead of 2h/m 
associated with the control cycle used in our calculation. 

It should be pointed out that according to our analy- 
sis in Sec. IIII Al the optical plug employed in the vortex 
pump simulation should render all K-quantum vortices 
dynamically stable at g = 250 up to k — 42, but here the 
splitting was observed already before reaching n = 22. 
The discrepancy is likely due to the fact that the stability 
analysis of Sec. IIII Al concerns a perfectly spin-polarized 
BEC in a pure harmonic trap. Thus, it can be strictly 
applied to the spinor BEC in vortex pumping only when 
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the external magnetic field consists of the strong homo- 
geneous bias field, i.e., at the start and middle point of 
each pumping cycle. The additional instabilities related 
to the spinor nature of the BEC and the presence of the 
multipole magnetic field have not been taken into ac- 
count, and hence the stability analysis of Sec. IIIIAl only 
provides lower limits of the stabilizing plug parameters 
for the vortex pump [56j . Nonetheless, the qualitative 
behavior A m - m oc cxp(aft) for some constant a is still 
anticipated in the spinor BEC. 
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